Helical structures from an isotropic homopolymer model 
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Abstract 

We present Monte Carlo simulation results for square-well homopolymers at a series of bond lengths. 
Although the model contains only isotropic pairwise interactions, under appropriate conditions this system 
shows spontaneous chiral symmetry breaking, where the chain exists in either a left- or a right-handed heli- 
cal structure. We investigate how this behavior depends upon the ratio between bond length and monomer 
radius. 
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Naturally occurring polymers, such as proteins (polypeptides) and nucleic acids (polynu- 
cleotides), are capable of adopting very specific single-chain conformations under certain con- 
ditions. It is believed that this structural specificity plays a crucial role in the biological function 
of a given molecule. Consequently, the essential mechanism of how the overall structure of a 
polymer arises from the basic interactions between its constituent monomers is of great funda- 
mental interest. A recurring motif in many diverse biopolymers is the helix. This structure has 
also been observed in synthetic polymers — for example, meta-phenylacetylene oligomers |Q], 
and the isotactic form of poly(methyl methacrylate) [2] in various solvents. The occurrence of 
helical structures in polymers with substantially different monomer chemistries suggests an un- 
derlying generic behavior, which has prompted significant efforts to produce "minimal" models 



that naturally display helical structures ||3|,|4|,|5|,|y]- An underlying assumption in the vast majority 
of these models is stated explicitly by Kemp and Chen [7]: "an isotropic potential interaction is 
not sufficient to produce helical ground states." In this letter, we present simulation results for a 
polymer model using only pairwise isotropic monomer interactions. We explicitly demonstrate 
that this model breaks chiral symmetry at low temperatures and show evidence that it forms stable 
helices. 

The general theoretical perspective 1^,0] has been to view helix formation as a discrete process. 
Each polymer element (monomer) can only take on one of two states, "helical" or "coil" (non- 
helical), and a free energy is assigned to each of these states. This approach has been fairly 
successful in describing the helix-coil transition in polymer systems. One shortcoming of this 
approach, however, is that it presupposes the existence of a particular helix, and, as a consequence, 
it is difficult to determine the conditions under which a helical conformation can naturally arise. 

Another line of research has been to apply molecular simulations techniques to examine the 
structural behavior of minimal polymer models, which has lead to interesting results. Zhou and 
coworkers discovered Jioiflli a freezing transition for linear homopolymers composed of tangent- 
hard spheres that interact with each other through an attractive square-well potential. During this 
discontinuous transition, the collapsed "molten" globule state freezes into a highly symmetric 
crystalline-like state. While this transition does result in structural specificity, the observed or- 
dered state does not resemble any of the structural motifs seen in biopolymers. Noguchi and 
Yoshikawa J12I extended the model of Zhou and coworkers by adding a bond bending poten- 
tial, which introduced stiffness to the polymer model. Upon increasing the stiffness of the chain, 
toroidal and rod-like phases were observed — analogous to structures exhibited by isolated DNA 
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chains JT3L IT4I1 . However, they did not observe the formation of helices in their simulation studies. 
Kemp and Chen Q] have developed a simple homopolymer model that displays a helix-coil 
transition; however, the monomer-monomer interactions in this model are asymmetric. 

For thick strings confined within a fixed geometry (e.g., inside a cube), a helix can be the 
optimal packing configuration ifisl . These "close-packed" configurations may play a role in de- 
termining the structure of linear chains in a similar manner as the close-packed structures for hard 
spheres play for the solid phase of spherically symmetric molecules. Structures other than the 
helix were shown to optimize the packing under appropriate conditions; the nature of the opti- 
mal structure depends on such parameters as the thickness and length of the string, as well as the 
geometry of the confinement. 

Although anisotropic, attractive interactions between monomers can induce a chain to form a 
helix, we will demonstrate that the delicate interplay between purely isotropic, attractive interac- 
tions and packing (due to excluded volume interactions) of monomer segments can also induce 
a helix. A homopolymer chain in poor solvent conditions will tend to collapse into a compact 
conformation in order to maximize self-interactions. Beyond a critical packing fraction, optimal 
configurations tend to maximize the entropy available to a system. Therefore, for situations where 
the helix is the optimal configuration for a string, we expect a transition to a helical structure for a 
sufficiently string-like polymer model. 

We consider polymers consisting of N linearly bonded hard spheres of diameter a. The distance 
between bonded spheres are restricted to lie between 0.9/ and 1.1/, where / is the bond length, 
which is not necessarily equal to a. Spheres that are not directly bonded interact via a square-well 
potential u(r), given by: 



u{r) 



00 for r < a 

— e for a < r < Act (1) 
for Act < r 



where r is the distance between the centers of the spheres, Act is the width of the attractive square- 
well (A = 1.5 in this work), and e is the magnitude of the interaction, which sets the energy scale 
of the model. For / = a, this model is precisely the same as that studied by Zhou and coworkers 

HQ. 

We have carried out Monte Carlo simulations to study the behavior of unconfined, isolated 
polymers with N = 20 (20-mers) and values of ct/7 between 1.0 and 1.9. Parallel tempering iflfil 
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was employed, with 16 replica systems at temperatures between /c^T/e = 0.2 and ksT/e = 7.0 
(where ks is the Boltzmann constant), placed such that the ratio between successive temperatures 

n 

was constant II17H . Within each replica system, the configuration of the polymer was modified 
using a combination of pivot moves 111811 . kink-jump moves [19], and continuum configurational 
bias regrowth moves H20I1 . Bond length fluctuations were implemented by performing standard 
particle displacement moves [21] on individual monomer spheres. After 400 attempted moves in 
each box, replica exchange moves were attempted. Samples were taken every 80 attempted swap 
moves. The simulation consisted of a total of 8 x 10 5 attempted swap moves. Smooth interpolation 
of data between the replica temperatures has been performed using multiple histogram reweighting 



[22] 



The structure of the polymer conformations has been characterized using the local torsion T; 
and curvature Ki [12311 around the bond vector joining monomers index i and i + 1 (i E [1,N]). 
Curvature is defined so that is always positive; the sign of the torsion around a bond indicates 
whether it has a left- or right-handed "twist". These quantities are calculated using a central 



difference scheme Il24h : as such, is defined on i E [3, N — 2], and Ki on i E [2, N — 1). It is 



a necessary and sufficient condition for a perfect helical structure to have Tj/ Ki independent of i 



1 25], taking the value p/2rrr h , where p is the pitch of the helix and r h is the radius. This quantity 



will be quoted as an average along the polymer for a given conformation, denoted 7, and given by: 

N-2 

7 = / ( N - 4 ) (2) 

This conformation average should not be confused with the usual thermal (ensemble) average. As 
an indicator of "helicity", we also measure the conformational root variance for this average, s 7 . 
This quantity should take consistently low values for a helical phase, since the ratio / Ki will be 
constant for a perfect helical configuration. 

The variation of the isochoric heat capacity C v with temperature T for a typical 20-mer chain 
(cr/Z = 1.6) is shown in Fig. [Ha). Cooling from high temperature, the heat capacity rises through 
a broad peak or shoulder as the radius of gyration (not shown) undergoes a large, though grad- 
ual, decrease. This is the "6-point" transition from extended coil to disordered globule, where 
the repulsive excluded volume interactions between the monomers are balanced by the attractive 
interactions. This transition is analogous to the vapor-liquid transition in bulk fluids. A snapshot 
configuration from the globule phase is shown in Fig. Eta). 
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Figure 1: (a) Variation with temperature of the isochoric heat capacity C v /k,B and mean conformational 
root variance, (s 7 ) for a homopolymer composed of 20 monomers, a/l = 1.6; jackknife estimated error 
is never more than 2% of the values shown, (inset) Mean end-to-end distance (r e ) for the same system; 
jackknife estimated error is never more than 2% of the values shown. Dotted drop lines indicate the tem- 
peratures of maxima in C v . (b) Probability density histograms P(*y) for the ratio of configuration average 
torsion and curvature for the same system, temperatures as in legend. 

Upon further cooling, the heat capacity passes through another peak at kgT/e = 0.58. We 
attribute this peak to the adoption of a compact, ordered structure by the polymer and denote 
this temperature Tf. Snapshot configurations taken from the simulations below this temperature 
appear to be helical in character. An example is given in Fig. [21(b). Since snapshot configurations 
are not a reliable indicator of structure, we have characterized the observed structures using the 
conformational mean torsion-curvature ratio 7, and its conformational root variance s 7 . 

Probability density plots P{^) for this system are shown in Fig.[T](b). At high temperature, 
these have an approximately Gaussian form, centered on zero. As the system is cooled below Tg, 
P{"i) gradually becomes double peaked, with the two peaks symmetrically placed around zero and 
separated by a shallow trough. The free energy barrier between states with 7 > (right-handed) 
and 7 < (left-handed) is relatively low, and the polymers can freely exchange between them. 
There is no clear thermodynamic indicator of the temperature below which this characteristic 
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(a) (b) (c) (d) 



Figure 2: Snapshot homopolymer configurations (a) 20-mer with a/l = 1.6 and ksT/e = 0.75 (globule 
phase), (b) 20-mer with a/l = 1.6 and ksT/e = 0.48 (helix 1 phase), (c) 20-mer with a/l = 1.6 and 
k B T/e = 0.2, (helix 2 phase) and (d) 32-mer with a/l = 1.9 and k B T/e = 0.2. 

appears. 

Below Tf, the character of P("f) changes once again. As the temperature passes below Tf, the 
two symmetrically placed peaks become separated by a region of very low probability. This marks 
the breaking of chiral symmetry, with the system adopting either a counterclockwise (7 < 0) 
or a clockwise (7 > 0) helix. These two states are separated by a large free energy barrier, 
and, consequently, the polymers are not expected to exchange between the counterclockwise and 
clockwise helix conformations. It should be noted that sampling of both peaks occurs mainly 
because of the use of parallel tempering in the simulations. That the two peaks are sampled with 
approximately equal weight is evidence that the simulations are ergodic and well equilibrated. 
Passing down in temperature through Tf, a sigmoidal decrease in (s 7 ) to a low, approximately 
temperature-independent value is also observed, showing that the t/k values for individual bonds 
are becoming strongly correlated, indicating a transition to a helical state. 

It can be seen from Fig.[T](a) that there is a further low-temperature peak in C v at ksT/e = 0.3. 
Upon cooling across this temperature, P (7) remains double peaked and shows a sudden shift 
inward, toward zero. Snapshot configurations remain helical, with an example shown in Fig.|2](c). 
The inset to Fig.[l](a) shows that the mean end-to-end difference for the polymer decreases upon 
cooling across this temperature. This tallies with the observed decrease in 7, indicating a decrease 
in pitch and/or increase in radius for a helical configuration. 

The general, qualitative behavior shown in Fig. |T| is also observed for polymers with a/l = 
1.475 to a/l = 1.65. For values a/l > 1.7, the behavior is as described above, but with only 
a single low-temperature peak observed in the isochoric heat capacity. For values a/l < 1.475, 
chiral symmetry breaking is not observed; the observed histograms for P(y) remain unimodal and 
centered on zero across the studied temperature range. These observations are summarized in the 
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Figure 3: Schematic phase diagram. Points represent observed isochoric heat capacity peaks; the statistical 
uncertainty in the positions is smaller than symbol size. Lines serve only as guides to the eye; see text for 
discussion. Text indicates the structure observed in that region. 

schematic phase diagram presented Fig.|3J 

We have shown that stiffness of the polymer induced by the bulkiness of the monomers plays 
an important role in the development of structure at low temperatures. For values of a /I > 1.475, 
the equilibrium low temperature structure appears to be helical (see Fig. El); further evidence for 
this is provided by the consistently low values for the quantity (s 7 ) under these conditions, indi- 
cating an approximately constant value for 7i//Cj across individual conformations. The torsional 
symmetry breaking associated with helix structure formation must be associated with a thermo- 
dynamic phase transition for an infinite system; no attempt is made to establish the order of this 
transition, as the small size of the system studied would make any such attribution unreliable. We 
have identified these transitions by the associated increase in energy fluctuations (maxima in heat 
capacity). Across the temperature range studied, no ordered low temperature structures are ob- 
served for all < 1.475. For a /I = 1, the globule should freeze into an ordered crystalline phase 

On 
lllll at sufficiently low temperature. For other values of a /I, however, the low temperature 

structure is still unclear. 

For values 1.475 < a /I < 1.66, the helix-globule transition line is reentrant. With increas- 
ing a /I, the helix initially gains, then rapidly loses stability with respect to the globule. At low 
temperature, there is a line of heat capacity maxima which coincides with a sudden drop in the 
torsion-curvature ratio and decrease in mean end-to-end distance for the observed helices. We in- 
terpret it as indicating a transition to a less extended helical structure at low temperature, since t/k 
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gives the pitch-radius ratio for a helix. This line of maxima appears to join with the helix-globule 
transition line for a/l > 1.675. No further heat capacity maxima or discontinuities in torsion are 
observed between this line and the lowest temperature studied. Since the high temperature helix is 
not observed for a/l > 1.675, and there must exist a transition line between a symmetry breaking 
and a non-symmetry breaking state, we can therefore infer the existence of a helix-helix-globule 
triple point in the range 1.66 < a/l < 1.675. The less extended helical structure is the low- 
temperature equilibrium configuration observed for all a/l > 1.475, and appears to increase in 
stability monotonically with a/l. Extrapolation of the helix-globule transition line and the helix- 
helix transition line to low temperatures for a/l < 1.475 suggests that a further, low temperature 
helix-helix-globule triple point may exist. 

The precise reasons for the occurrence of two helical structures separated by an apparent transi- 
tion line, and for the sudden decrease in stability of the high temperature helix are, as yet, unclear. 
Results using an anisotropic model JsIlL constructed to yield a specific helix, suggest a two-step 
mechanism for helix formation, where the initial globule-helix transition results in a structure with 
disordered end segments, followed by a transition to a fully-helical structure. Another speculative 
possibility is that allowed helical states are discretized due to a preference for monomers to sit 
at interstitial sites in neighbouring layers of the helix, leading to a succession of helical phases 
separated by transitions rather than a continuum of possible torsion values. Further study will be 
necessary to understand these behaviors. 

Another factor playing an important role is the overall length of the polymer. We have studied 
the square-well chain model for different lengths N, and preliminary results suggest that, though 
the system breaks chiral symmetry at low temperature for all chain lengths studied, helices are 
only stable inside a window of values for N. An example low-temperature snapshot configuration 
for a longer polymer is shown in Fig.|2](d); a detailed examination of the influence of chain length 
on structure will be the subject of a further communication. 

In summary, we have demonstrated spontaneous chiral symmetry breaking using a model which 
is achiral (no dependence upon bond torsions) and isotropic (no dependence upon bond angles). 
This spontaneous symmetry breaking reflects the formation of helical structures at low tempera- 
ture, providing support for the notion of helices as "optimal configurations" for string-like objects 
I15I1 . which as such do not need to be "engineered in" to a model or molecule through specific 
directional interactions. We do not seek to suggest that directional interactions do not play an 
important role in the formation and stabilisation of helical structure; however, the purely steric 
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effect which we have demonstrated may also play a crucial role which we have not yet been aware 
of. The observed phase behavior for the model is rich, including an apparent helix-helix transi- 
tion, with a reentrant helix-globule coexistence curve for one helix leading to a helix-globule-helix 
triple point, and offers many possibilities for further investigation. 
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